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We study the finite temperature behavior of the CPT-even pure-photon sector of the Stan¬ 
dard Model Extension, which is defined by the standard Maxwell Lagrangian plus the term 
{kF) l u/cipF IJ ' 1 'F a P. The Hamiltonian analysis is performed, from which the degrees of freedom and 
constraints of the theory are derived. We have explicitly calculated the partition function for an ar¬ 
bitrary configuration of the (kp)^^ coefficients, to second order, and we have used it to obtain the 
thermodynamic properties of the modified photon sector. We find the correction to the frequency 
dependence in Planck’s radiation law, and we identify that the total energy density is adjusted, rela¬ 
tive to the standard scenario, by a global proportionality constant containing the Lorentz-violating 
contributions. Nevertheless, the equation of state is not affected by these modifications. 
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I. INTRODUCTION 

One of the cornerstones of quantum field theories and general relativity is Lorentz invariance, which is assumed to 
be an exact symmetry. There is a lot of remarkably strong experimental support for this idea flj, with no violation 
detected. Nevertheless, certain quantum-gravity theories possess mechanisms that can lead to Lorentz violation Q, 
in which effects arising from modifications of space-time are expected to appear at distances of the order of the 
Planck length. This has attracted considerable attention in recent years both from the experimental and theoretical 
perspectives, given that, in principle, this would allow a better understanding about the space-time structure. The 
Standard Model Extension (SME) performed by Kostelecky et al. Q is a framework where Lorentz violation is 
motivated via spontaneous symmetry breaking (SSB), in which non-dynamical tensor fields are now added to General 
Relativity and the Standard Model and whose fixed directions induce the corresponding symmetry breaking in a given 
reference frame. These tensor fields are assumed to arise from non zero vacuum expectation values of some basic fields 
belonging to a more fundamental model, such as string theory Q. Considering the SME as a framework motivated 
from a SSB, the Goldstone theorem ensures that massless particles will emerge Q. Bjorken proposed that the photon 
can be a Goldstone mode associated with the SSB of Lorentz invariance [g]. However, since the pure-photon sector 
of the SME is 17(1) gauge invariant, all particles are massless, and an alternative interpretation of the Goldstone 
theorem is required [7|. In fact, in the spontaneously broken space-time symmetries case, the counting of massless 
modes has to be done carefully j§], and it may happen that one to six Goldstone modes appear, each corresponding 
to one of the six Lorentz generators. The properties of Goldstone modes are, in general, model dependent and the 
knowledge of the fundamental theory is required to do a complete description. In the case of the SME some general 
conclusions about Goldstone modes have been obtained in Q. 

The pure-photon sector of the SME includes the usual Maxwell term plus the presence of the CPT- 
odd term {\{kAF) K ^K sometimes called the Carroll-Field-Jackiw term 0 , and the CPT-even term 
—\{} i F)nvapF 111 'F a P . Both terms have been extensively studied in the literature EL and experimental constraints 
exist for them jTsj. The search for new effects arising from these Lorentz violating terms, and an improvement of 
the bounds for the magnitude of these coefficients constitute two of the main lines of study. The study of the cosmic 
microwave background (CMB) offers an opportunity to test the pure-photon sector of the SME at finite temperature 
EL since the propagation of light would be affected^in the form of non-standard dispersion relations, polarization, 
birefringence properties, among other effects [1, 11ll- fl3)| . As it may be expected, thermodynamic properties and the 
spectral distribution can be modified as well. In Ref. [15[ , the partition function in the functional integral formalism 
was calculated in order to study the finite temperature behavior of the Carroll-Field-Jackiw term, for the case of a 
purely spacelike background. In Ref. EJ the study was extend to the CPT-even term for particular configurations of 
the coefficients (/cf)/x va/P to simplify the calculations. In both cases, Lorentz violation corrections to the blackbody 
radiation and anisotropy in the angular distribution for the energy density were found. It remains unclear, however, 
if there is no information loss from the consideration of only a few particular configurations. The goal of this paper 


2 


is to study the finite temperature properties of the CPT-even pure-photon sector of the SME in the most general 
case, following a scheme similar to that employed in [l5|, [n|. The outline of this work is the following. In Section ITT] 
the CPT-even pure-photon sector of the SME Lagrangian is introduced and some properties are reviewed. Following 
Dirac’s scheme for constrained systems, the Hamiltonian analysis is performed, and the canonical quantization is 
carried out. In Section uni the partition function is evaluated in the functional formalism and some thermodynamical 
properties are derived, a substantial improvement over the previous calculations is reported considering an arbitrary 
configuration for the ( k F )^ va ^ coefficients. We compare the result with a classical thermodynamic approach. Our 
summary and conclusions are contained in Sec. hyi 


II. THE MODEL 

We focus on the pure-photon sector, and particularly on the CPT-even violating terms within the minimal SME. 
The Lagrangian density is given by 

L = —F^F^ - \{k F )^ a pF^F a P , (1) 

where the first term corresponds to the standard electrodynamics, being = d ll A v —d^A^ l the electromagnetic stress 
tensor. The second term introduces Lorentz-CPT symmetry breakdown, which is controlled by the non-dynamical 
spacetime-constant and dimensionless quantities {kp)^^-, these coefficients have the symmetries of the Riemann 
tensor and a vanishing double trace, which imply a total of 19 independent components. The tensor (k F ) / j, va p is 
alternatively parametrized in terms of four 3x3 matrices 0 , KDE, khb,kdb,khe, defined by 

(KD E y k = - 2 (k F ) ojok , ( K HB ) 3k = l e jp<i e kim (k F r im , {K DB y k = -{K HE ) kj = ( 2 ) 

The matrices kde and khb contain together 11 independent components, while kdb and khe possess together 8 
components, which encompass the 19 independent elements of the tensor {k F )^ An alternative parametrization, 
which allows easier experimental constraints, consists of writing ( k F ) fJjV , a p in terms of four traceless matrices and one 
trace element (llj . 

{k 0 +y k = -(kdb + kheY 1 * , (k 0 -) jk = -^DB-K HE y k , (3) 

(k e+ ) jk = -(kde + khb)^, (k e ~y k = ^ KDE ~ K HBY k - -^S^Tt^kde), «tr = -Tr (kde)- ( 4 ) 

All parity-even coefficients are contained in k e+ ,k e - and k tr , while all parity-odd coefficients are contained in k a+ 
and k 0 _. The matrix k a+ is antisymmetric while the remaining matrices are symmetric. In Section IlIIl we will use 
this second parametrization to express our main results. 

As we previously mentioned, the coefficients can be motivated by spontaneous breaking of Lorentz sym¬ 

metry (4], avoiding the issues of incompatibility in General Relativity present when an explicit Lorentz symmetry 
violation is introduced 0. The transformation A^ + <9 M A leaves the Lagrangian ([1]) invariant, and therefore, 

the gauge symmetry U (1) is preserved as in the Maxwell theory. We use the convention, Greek indices p,, v = 0,1, 2, 3, 
Latin indices i,j = 1,2,3 and the metric rj^ = (1, — 1, — 1, — 1). The Euler-Lagrange equations arising from the 
Lagrangian Q correspond to 

d v F^ + (k F )^ a pd u F a P =0. (5) 


The remaining Maxwell equations 

d,F pv = l -e^ ali dii F Qp = 0 , 


( 6 ) 


continue to hold as a consequence of defining F through the potential A M . As we previously stated, the propagation 
of light is modified due the presence of the vap coefficients, in this case the dispersion relation for the Lagrangian 
© is given by jll| 



OL 


^ = ^(ka0) 2 -p 2 , 


P°± = (1 + p±a)\p\, 


( 7 ) 
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to lowest order in (k F )^, va f3, with 


~k ap = (k F ) a ^ v p^, r = |j • (8) 

Let us now consider the analysis of constraints a la Dirac flsl , and the canonical quantization of the model. This 
will allow us to determine the number of degrees of freedom and establish some differences between the standard 
electrodynamics and the CPT-even pure-photon sector of the SME. The canonically conjugated momenta are given 
by 


7r‘ = -y = F M + (k F y 0a ^F a p (9) 

oAi 

= f>\5 i i + 2 (k F y° j0 ) + ( k F y olm F lm 

= M i j F j0 + N i , 

where we have defined M l j = 6 l j + 2(k F ) M jo and N l = ( k F ) l0lTn Fi m . The canonical momentum associated to Aq is 
null, 7T° = 0, as it is in standard electrodynamics. Approximating all quantities to first order in the Lorentz-violating 
coefficients, we find that the matrix M % j has the inverse B x j = (M -1 ) 1 j = d l j —2 (k F ) M jo, which allows us to rewrite 
m as F k0 = B k iir 1 — N k . Using the above it is straightforward to obtain the canonical Hamiltonian density 

H c = -i B 1 j 7r^7r j + -F l3 F^ - {k F ) 0ilm Fi F lm + l(k F r lm F tJ F lm - A 0 d k n k , (10) 

where we have carried out one integration by parts and omitted boundary terms. The non-zero Poisson brackets (PB) 
are given by 


(4u( x > *)> ^(y. t)}p = v 6 3 (x - y). (11) 

In what follows we will assume that all PB are calculated at equal times and we omit the label t. We employ Dirac’s 
method to construct the canonical theory due to the fact that the primary constraint 

~ 0 , ( 12 ) 

is present (here the symbol ~ denotes the weak equality). The extended Hamiltonian density is defined as 

H e = H c + A0i, (13) 

where A is an arbitrary function. The evolution condition of the primary constraint (El). 


0i(x) = (0i(x), J d 3 y H E (y)} P ~ 0 , 


(14) 


leads to Gauss’ law, 


02 = di-K 1 ~ 0. (15) 

It is not difficult to prove that El and El are the only constraints present in the model, and that they con¬ 
stitute a first class set ({0i, 02} = 0). As in standard electrodynamics, the model possesses two degrees of free¬ 
dom (DOF), following Dirac’s scheme, DOF = |[variables in the phase space — second class constrictions — 2 x 
first class constrictions] = ^[8 — 0 — 2 x 2] = 2. If we write (fT5l) in terms of the field strength using the defini¬ 
tion of the canonical momenta ©, we obtain 

di^ = &F 10 + {kp^diFcp = 0. (16) 


We recognize the last equation as the Lagrangian equation © for the /i = 0 component; however, in the Hamiltonian 
formalism it is a constraint, not an equation of motion. Gauss’ law El is different from its standard electrodynamics 
form, being this even more evident when it is rewritten in terms of the electric and magnetic fields instead of the 
canonical momenta (7P = F t0 + (k F y 0a/3 F a p). In order to construct a quantum theory via canonical quantization 
({A, B} —>• A-[A,B]), we must remove the extra degrees of freedom. This means that we have to impose as many 
suitable gauge constraints “by hand” as there are first class constraints; these gauge constraints have to be admissible 
and convert the first class constraints into second class constraints, and then we can introduce the Dirac brackets to 
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perform the correct quantization. We choose the Coulomb gauge (4>i = diA 1 ~ 0) plus $ 2 = A 0 ~ 0 to fix the gauge. 
The Dirac brackets 


{A(x),B(y)} D = {A(x),B(y)} P - 


d 3 ud 3 v{A(x),Xi(u)} P (Q ^ {Xj(v), B(y)} P , 


(17) 


where \i is one of the constraints </> 2 , $i, $2) and Q y (x, y) = {%j(x), Xj(y)}P: 

( 


Q l ° (x, y) = 


are given by 


( 0 0 -10 
0 0 0 V 2 

10 0 0 

Vo -v 2 0 0 


(x - y); (Q~ 3 ) v (x,y) = 


0 
0 

—<5 3 (x — y) 

0 


0 

0 

0 

1 

47r|x—y | 


<5 3 (x-y) 0 ^ 

*' 4 tt|x— y| 

0 0 

0 0 ) 


M*(x,i),A J (y,i)} D = 0, 
{n l {x,t) 1 n J (y 1 t)}D = 0, 

{Ai(x,t),ir 3 (y,t)}D = Si 3 S(x - y) + did 3 


47t|x — y| 


= ^L^x-y), 


(18) 

(19) 


and they have the same form as in standard electrodynamics when the canonical momenta are used. However, 
rewriting in terms of the electric and magnetic fields, the difference is manifest, 


{A 4 (x, t), A 3 (y, t)} D = 0, (20) 

{E i (x,t),E 3 (y 1 t)}D = 2 [(k F y oh d x i - ( k F ) 3 ° h d y i}6 3 {x - y), 

{A a (x, f), E 3 (y, t)} D = 3 (x - y) + 2 (k F ) 3 ° s o^L s (x - y). 


Once the Dirac brackets have been included, the dynamics of the theory will be generated by the Hamiltonian (fTUl) 
without the term A 0 diTr z , which is proportional to </>2 and has already been fixed and included in the process of 
introducing the Dirac brackets. The canonical quantization can be carried out using the aforementioned Hamiltonian 
m and the brackets given by m- Being 17(1) the group of symmetry of the theory, there are other possibilities 
to fix the gauge, as the Lorentz gauge (9 M A M =0), which is manifestly covariant; however, it is not possible to handle 
such a gauge with the quantization formalism that we employ here. This is due to the fact that the Lorentz gauge 
involves the time derivative of Aq. There are well known formalisms which are capable of dealing with such relativistic 
constraints, among which are BRST quantization [l9j, the Fadeev-Popov method [2(j] into the path integral (2l[ or 
within the Hamiltonian formalism, one has [Hj]. Nevertheless, these approaches are not within of the scope of the 
present work. 


Rewriting the Hamiltonian (flUl) in terms of the electric and magnetic fields we find 

H = i(E 2 + B 2 ) - (, k F ) 0j0k E 3 E k + i( k F ) 3klrn e 3kp e lrnq B p B q , (21) 

where E 2 = EiEi, Bi = —±e l3k F 3k , and therefore B 2 = ^F^F 13 . The same result has been found in [3j following a 
different line of thought, where it was shown that if (kp )^ l/a ^3 is small, the last quantity (1211) is nonnegative. This 
is due to the fact that the Hamiltonian m can be viewed as the bilinear form x T Mx with x T = (E,B). It can be 
readily shown that, upon diagonalization, the matrix M has entries 7 — 0(kp) > 0 for both the electric and magnetic 
field contributions. 


III. PARTITION FUNCTION AND THERMODYNAMICS 


We derive now some of the thermodynamic properties of the Lagrangian 0. Our main goal is to obtain the finite 
temperature energy density of the electromagnetic field. Following the quantum field theory scheme, we calculate the 
partition function. In the previous section we adopted the Coulomb gauge; hereafter we will switch to a covariant 
gauge. The simplest way to obtain the partition function is trough the Faddeev-Popov method [20j, which is equivalent 
to the introduction of constraints as done in Sec. El both methods allow us to work with the effective degrees of 
freedom. Choosing the Lorentz gauge we can write the partition function as 


Z = 


J[dA ^ det(-d 2 ) exp 




( 22 ) 
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where det(— d 2 ) is the Faddeev-Popov determinant and we have switched to an imaginary time variable t = it. The 
effective Lagrangian is given by 

L eS = -\f^ v ~-^k F )^ a pF^F^ - y p {d^f . (23) 

Upon substitution of the finite temperature replacements t —> —*r, Aq —>■ iA T and ( k F ) 0llva —>• i(k F ) T,lva (simi¬ 
lar convention to other indices in (k F ) tll,0 ‘P), the effective Lagrangian can be written in Euclidean notation, with 
ii,v,a,p = r, 1,2,3, 


Leff — 




5^d 2 - ( 1 - - )d tl d v + 2 {k F )p vaF dpd 0 


A„ 


= ^ 

1 T 

= -A t DA. 
2 


$ [IV d A 2 (k F ')pva^dpdQ 


(P^l) 


(24) 


In the first line we have chosen the Feynman gauge (p = 1) and D^ u = S^d 2 + 2{k F )p l , all dpd a . The field admits a 
Fourier expansion: 


A m (t, x) 



n ,p 


e^+X'P } i>,p), 


(25) 


where ui n = are the Matsubara frequencies, and the field A m (t, x) satisfies the constraints of periodicity A^fi, x) = 
A m (0,x) for all x. The normalization in (l25l) is chosen so that each Fourier amplitude is dimensionless. If we use a 
ghost field C to write 


det(— d 2 ) = J [dC][dG]exp yj dr j d 3 x(d l _ l C)(d ,J 'C) J, 
then we can calculate the partition function in frequency-momentum space as 

In Z = Trln[/3 2 (w 2 + p 2 )] - iln[Det(D)], 


where now 


( ^n IP t A rr A tx 


D = 


, 2 ^p 2 

A 
A. 
A. 


TX 

ry 


< + p‘ 
A x 
A, 


A 

A x 


ry 


*-xx l *-xy 

wl + P 2 -+ 
A yz 


TZ 

XZ 


w n + P“ + A- 


and A^ v = 2(k F )[ lav pp a pp, ( p T = w n ). Calculation of the determinant to second order in k F gives 
det(D) = \{^[{u 2 n + p 2 ) 4 ] fl + Tr(A) + i(Tr(A)) 2 - ^Tr(A 2 )) , 

n,p ^ ' 

where we have defined A ^ = 2{k F )^a^pPaPp /(w 2 + p 2 ) and the relation 


Det(l + M) = 1 + Tr(M) + ^(Tr(M)) 2 - ^Tr(M 2 ) + 0(M 3 ), 


(26) 


(27) 


(28) 


(29) 


(30) 


has been employed. Therefore, the total partition function becomes 


In Z 


- £ ln[/3> 2 + p 2 )] \Jl + Tr(A) + i(Tr(A )) 2 - ^Tr(A 2 ) + • • ■ 

n,p n, p ' 

In Z 0 + In Z lv ■ 


(31) 
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We recognize the first term as the usual result for the Maxwell theory, which corresponds to massless bosons with 
two spin degrees of freedom in thermal equilibrium; in other words, blackbody radiation. All modifications to the 
standard case due to Lorentz violation come from to the second term in (13111 . which we can evaluate as follows 

In Zlv = ln(l + ^(A) + i(W(A)) 2 -iw(A 2 )') » -± £ (w(A) - Er(A 2 )) =-^(Z IVi +Z I% ). (32) 

n,p ' ' n. p ' ' n, p 

Here we have defined Z F y x and Z F y 2 as the Lorentz violation contributions to first and second order, respectively. 
We begin calculating the first order contributions as follows: 


= -^E Tr ( A ) 


n.p 


n.p 


= -E 


{ L k F ')a F oti'PiJ.Pu 


+ P 2 ) 


27t n. 

[Po =u n = —) 


-Y _ 1 

Z—l 2 4- 


n,p 


(^ + p 2 ) 


~Y _- 

C /,,2 


n, p 


(^ + p 2 ) 


{k F ) OLT olt + 2(/cf) OLlOLT MnPi + {^F^aiajPiPj 

(k F ) 

OLTOLT w 2 + (. k F ) aiaj piPj 


(33) 


where in the third line we employed the sum 'F'^L-oo n' 2 +a 2 = 0- Adding and subtracting the term (k F ) aTaT p 2 inside 
the brackets in (1331) . we arrive to the equivalent expression 


-J2 Zlv i = -E 


{k F ) c 


— (k F ) c 


( k F ) aiajPiPj \P 


21 


(w 2 +P 2 ) 


(34) 


with pi = Pi/\p\- Making use of the identity 

“ 1 


E 


= — coth 
<^n + P" 2 p 


/3p\ _£ 

2 p 


1 + 


;/3p _ 1 


(35) 


we find that the contribution due to Lorentz violation to first order can be written as 


~Yz LVi = -v 


d 3 p 

(2f)3 


[~(k F )c 


(k F ) aiajpiP j]P 


2 p 


1 + 


iPp - 1 


(36) 


where a temperature-independent divergent term has been dropped; it is well known that any quantity in finite 
temperature theory is defined after subtraction by its T = 0 counterpart, and so T -independent parts of the partition 
function, infinite or finite, are of no importance 23] . Taking the standard spherical coordinate system and |p| = w, 
eq. m becomes 

~ ^ ' Z LVl = — J du) dil [— (k F ) aTolT + (k F )aiajPiPj\ ^ ^ j ^1 f ^ > (37) 


where pi = sin(0) cos(0), p 2 = sin(0) sin(^>) and P 3 = cos (6). Performing the integral over solid angle we find 


- E Zlv i = ~ 4?r 


n,p 


v 1 

771 77 [ ( k F )arar T 7j(.k F ^) a iai] 
[2 n)° o 


did 


1 3u.r 


eP u - 1 


(38) 


We now make use of the vanishing of the double trace of {k F )^ va p condition, which in Euclidean space implies 

(k F Y v = 2 {k F )° l o i + ( k F ) lj ij = —2(k F )om + ( k F )mj ———* 2 (k F ) TiTi + ( k F )ijij = 0 . (39) 

Using the above, (1551) becomes 

r du (Y-)(i + -Z-Z] = 2 AL(M™ = (40) 


\ ^ ^ 167T V n ^ 

- £ Zl Vi = ~ 5 ~ 777777 ( k F )c 


n,p 


3 (2tt) 3 


.fiui _ 1 


45/3 3 


45/3 3 tr ’ 
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where we have neglected vacuum contributions to perform the integral. The above implies that the modification to the 
energy density due Lorentz violation to first order will have the same dependence in the temperature as the standard 
theory, U ~ T 4 . The second order contribution can be readily evaluated in a similar way, 




<4 . , P 2 ^n . , P 4 

1 (^+p2 ) 2 + 2 (w 2 +p2)2 + 3 (w 2 +p2)2 > 


(41) 


where the momentum-dependent coefficients ^ 1 , 2,3 are defined as 

A\ — jiTVT {_^F^) firur •> 

A2 — lltvt (J^f')\ i]^F^)^iTvi H - ^f') K^f') firisj H“ ( &F ) (ij vt]^P iPj j 

A3 = {^f) fiiuj {^f) filiymPiPjPlPm • 

Making use of the identities 


(42) 


E 


n =—00 


44+ P 2 ) 2 8p 2 


^ crftlPl' 9 


41 p | a 


1 + 


e /3|p| — 1 / ’ 


(43) 


OO 9 

E —^ 

n= — 00 


44+ p 2 


= -?csch 2 


IPl/3 


4|p| \ e^l p l — 1 J 


(44) 


E 


= E 1 - 


2^P 2 


_(^+P 2 ) 2 4^ V 44+PT 44+P 2 )4 

we now evaluate the sums over n. Employing the same spherical coordinate system as in (EH) we find 


-E^% = v 


= V 
= V 


d 3 p 

(2^)3 

d£lduj 

du 
2 n 2 


{Ai - A 2 + yl 3 )EA csc h 2 ) + (-371! + A 2 + 7l 3 ) 


, . . . . w 4 /? 2 

{Ax - A 2 + A 3 ) csch 


2 

21 uP 


W( 


1 + 


+ (-3,41 + A 2 + zl 3 ) 


4 \ e^lPi - 1 

w 3 /3 /. 2 


1 


e /3w _ 1 

2 


~ ~ - w 4 /3 2 2 / w/3\ . - - - w 3 /3 

(Ai - A 2 + A 3 )—csch f — j + (-3^i + A 2 + A 3 )^—I 1 , e/3aj _ 1 


(45) 


(46) 


where now the momentum-independent coefficients A correspond to 

Ai = Ai = -Ti-(k 2 de ), 

9 / 

= | (M /LTVT 4.F ) fiiisi 4“ (4f) flTviiJ^F^) fiTisi 4“ ) fj,Tvi {kp ) fiiuT 


Tr (kde) ~ ■ khb) — 3Tr(ft BB • khe) — Tt(kdeY 


A 3 — _^ r ( (^f) fiivi(kp)lijvj 4“ (kF')/j,ii'j(kF')iiivj 4“ (&f) fiivj {^f)x 


1 

30 


Tr (kde • fde) — Tt(kdb • kdb) — 4Tr(/«£> b • khe) 4- -Tr(K ffB • khb) 4- Tr(« BB ) 2 


(47) 

(48) 

(49) 

(50) 


The integral over w can be calculated neglecting vacuum contributions, which arise from the second term in (1461) and 
have the same form as they do to first order in kp- The result is given by 

-J 2 z lv 2 = v {i A ^- A2+A ^4^i^) + i ~ 3Al+A2+A3) ^{rk i )) 


(51) 






































= K 


45/3 3 


1/, 


(52) 


where we have defined 


I< =-(Ai ^ 3 A 2 + 5 A 3 ) 


1 

16 


8Tr(«D£;) 2 — 2Tr (kdb ■ kdb) + 10Tr(/t db ■ khe) — Tr (kde • kde) 
+ 6Ti(k de ■ khb) + 7Tt(k H b ■ k H b) 

2Tt(k 0 + • k 0 +) - 3Tr(« 0 _ ■ k 0 ~) - Tr(k 0+ ■ «o-) + 3Tr(re e+ ■ k e +) 

- 4Tr (k e+ ■ k e ~) - 4(K tr )Tr(K 0 _) + 18(K tr ) 2 


(53) 


The partition function of the standard Maxwell theory is well known j24j : neglecting the vacuum contributions it is 
given by 


In Z n = -2V 


d 3 p . a , ,, TT“ 


(2tt) 


rln(l - e~ pu ) = 


45/3 3 


rV. 


From these results we finally obtain the total partition function 


In Z — In Zo + In Z E y — (1 + 2 (kp) aTaT + K) „ V. 

45 


(54) 


(55) 


In order to compare this result with the literature we consider the particular case arising from the isotropic contribution 
of the parity-even sector, which corresponds to the limit Tr(re 0 _) = Tr(«; 0 _|_) = Tr(«; e+ ) = Tr(« e _) = 0, k tI ^ 0. In 
this case (l55l) reduces to 


111 Z — ( 1 + 3(«tr) + ^(Ktr)" 


45/3 3 


V. 


(56) 


To second order this is the same result reported in [If 


Starting from (@qd, dH and dSU), the energy density of the photon field can be calculated from the standard 
thermodynamic relations, 


1 din Z 


V dp 


(57) 


did 


2 (kjr)c 


1 id J 
7T 2 — 1 37T 2 

r °° did 


did 


fhd^e 


4 


+ (-^■1 ~ A 2 + 243 ) 
~~( — 3 j 4 i + A 2 + A 3 ) 


0 !67t 2 

r °° chd 


_eP u — 1 {eP u - l) 2 
^‘(/Jacoth(^)-2lc s ch 2 (^') 


/ 0 


4-7T 2 




4 


e0u _ 1 ( e /3w _ 1)2 


(58) 


This implies a modification to the Planck distribution, where now the frequency dependence of the energy density is 
given by 


i(td) = 


1 w 3 

7T 2 — 1 7T 


1/2 1 

— (/cp’) QraT + -(— 3 Ai + A2 + A3) 


id 


/3uj 4 e 


4„0u> 


+ 16?r 2 + ^ /duj4 ( coth 


pid 

T 


- 2 I 


csch 


e? u - 1 (e/ 3 ^ - l) 2 

pid 


(59) 


(see Fig. [TJ. The total energy density is obtained integrating the last equation or deriving directly (l55l) 

u(T) = (1 + 2 (k F ) aTaT + K)^T 4 . 


( 60 ) 
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The thermodynamic relations 


u =-w l * z - 




together with (1551) . imply that the equation of state remains unchanged, 

PV=l. 


(61) 


(62) 



Figure 1: Frequency dependence of the photon energy density for non-vanishing CPT-even violating terms, to first order in Icf- 


The functional scheme that we used to derive the photon energy density is not the only prescription available. It 
is possible to reproduce the previous results following a classical thermodynamic scheme, as we now show. Since this 
computation is not the main purpose of our study, we only calculate the energy density to first order in fcp. We follow 
a similar treatment as j25j . By adopting a thermal distribution for the photons and dispersion law 0, the spectral 
energy density per frequency (and per polarization) in the semiclassical phase space is given by 

u ( w ±)d w± = (2 * )3 ehu± ^ T _ 1 k 2 dk dn . (63) 

From 0, we have uj± = (1 + <5±)|p|, where S± = p ± cr. For the isotropic frequencies, one can immediately substitute 
f dfi —» 47r, however we now have an angular dependence for the frequency, u± = w±(0, </>). To first order in S± -C 1 
and using |fc| = ui, we find 


u(u!±) duj± 


1 ( Huj 3 \ J jrt , 1 ( hio 3 ti 2 w 4 e ^/fc B T ^ r j jr> 

(27r) 3 \e^/ k ^ T - 1 ) dUjd + (2tt) 3 \ e ^/ k B T - 1 k B T (eW^T -iy) d± U}d 


(64) 


If we sum over both modes under the assumption that each one contributes equally, we have that the Lorentz violation 
contribution to the energy density is given by 


[u(uj + )du} + ) L v + [u(uj-)duj-}LV 


l / hoj 3 

(27r) 3 y e ftu/k B T _ i 


h 2 w 4 e R “/ fcBT V 
k B T (e s “/^ T - l) 2 ) + ^ 


(65) 


Since 6+ + 5- =2 p = —k a a and 


k a 


— IT 


■> {kf) OLTOLT 


‘Zi^kF^jrjiPi ^F^aioijPiPj > 


( 66 ) 
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where a = r, 1, 2,3, upon integration of (1651) over the solid angle dfl we finally have 

[u(u})cIw]lv = [ u(ijj+)d(jj+]Lv + [■ u(w-)cLj-]lv (67) 

4tt /n ^ l n ^ J hco 3 4 1 h 2 w*e nu / kBT \ J 

~ "(2^ ((fcF)ar “ T " 3“ k B T (e s “/^ T - l ) 2 ) dUJ 

2 fn ^ J Huj 3 h 2 ^e^/ksT ^ 

“ _ 1 - k B T (eWfcBT _ !)2) du} ’ 


where in the second line we have used the condition of vanishing double trace in the imaginary time (1391) . Taking 
h = 1 and j3 = 1 /k B T, we find a result in agreement with (l59l) . 


IV. SUMMARY 

We have presented in this paper a study of the Hamiltonian formalism and the canonical quantization, via the Dirac 
formalism, of the CPT-even photon sector of the Lorentz-violating Standard Model Extension. We have found that 
the gauge freedom of standard electrodynamics is not lost in the CPT-even photon sector of the SME. Additionally, 
in analogy with the standard electromagnetic case, this model possesses first class constraints which restrict the size 
of phase space; however, the form of these constraints is different from the standard case when they are written in 
terms of the electric and magnetic fields. The gauge fixing procedure in the CPT-even photon sector of the SME does 
not significantly differ from the standard electrodynamics case. 

The partition function was explicitly calculated to first and second order in the Lorentz-violating parameters of the 
CPT-even photon sector of the SME, for an arbitrary configuration of the coefficients (kp) l j.iiaPt we have presented our 
results in terms of the defining (fcp)^^ coefficients as well as the two parametrizations © and J3J4]). Making use of 
the standard thermodynamic relations we found corrections to the blackbody radiation law. These include corrections 
with the same functional dependence on frequency, as well as corrections with a different functional dependence. The 
total temperature-dependent energy density receives a correction with the same temperature dependence as the 
standard case, but we found that the equation of state has the same form of the standard electrodynamic theory. 
This corresponds to a cosmological thermal history consistent with the standard scenario. The result was verified 
employing a classical thermodynamic scheme. We have obtained that all modifications arising from Lorentz violation 
to first order are proportional to (k B ) a 0 a 0 coefficients, there is no contribution to first order from the ( k B )oi a /3 sector 
to the partition function as well as the density energy. We have compared our general results with a particular case 
previously obtained in [lf| for the isotropic contribution of the parity-even sector, which is encoded into the coefficient 
(k B )iOio ~ ktr, finding a complete agreement between both works. Furthermore, following the scheme employed here, 
it is possible to obtain the contributions arising from the couplings of the distinct sectors of the full CPT-even photon 
sector of the SME, as explicitly shown in Eq. HMD- 
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